econ401 final project R Notebook

Libraries:

library(tidyverse)
library(readr)
library(rvest)
library(stats)
library(readxl)
library(dplyr)
library(stringr)
library(ggplot2)
library(ggthemes)
library(stringr)
library(data.table)
library(geojsonio)
library(leaflet)
library(rgdal)
library(haven)
library(stargazer)

Data:

Model 3:

First making a graph of annual GDP highs and lows – maybe a temporary proxy for recessions?

Annual enrollment graph:

Join enrollment data and gdp data to create linear model test:

Graph it?

College Proximity Question 5/3: (Reading in Ivy’s data)

cz <- read_dta('cz.dta')

Error in read_dta("cz.dta") : could not find function "read_dta"

Read in/create mobility data: (Trends in Mobility: Commuting Zone Intergenerational Mobility Estimates by Birth Cohort) https://opportunityinsights.org/data/?geographic_level=101&topic=0&paper_id=0#resource-listing

Read in geojson file:

Commuting zones on the map (cz.geojson@data) are in 1990s format. They need to be converted so our post-2000 data can be connected to the shapefiles: (https://www.ers.usda.gov/data-products/commuting-zones-and-labor-market-areas/)

Try to run some lms:

Ivy’s STATA code:

foreach x of varlist ncollege nfouryr nfouryrpriv npub nelite hascollege{
    
foreach y of varlist kfr_pooled_pooled_p1 kfr_pooled_pooled_p25 kfr_pooled_pooled_p50 kfr_pooled_pooled_p75 kfr_pooled_pooled_p100 {
    reg `y' `x' popdensity2010 med_hhinc2016, r 
    outreg2 using `x'_kfr, excel append ctitle(`y')
}

Variables of interest:

as.formula(paste0(yvar1, " ~ ", paste0(xvars1, collapse =  " + ")))

kfr_pooled_pooled_p1 ~ ncollege + nfouryr + nfouryrpriv + npub + 
    nelite + hascollege + popdensity2010 + med_hhinc2016

testing?

stargazer(lm.kfr_p1,
          type = "text",
          dep.var.labels = c("kfr_pooled_pooled_p1"))

length of NULL cannot be changedlength of NULL cannot be changedlength of NULL cannot be changedlength of NULL cannot be changedlength of NULL cannot be changed

===============================================
                        Dependent variable:    
                    ---------------------------
                                kfr            
-----------------------------------------------
ncollege                       0.003           
                              (0.004)          
                                               
nfouryr                       0.007*           
                              (0.004)          
                                               
nfouryrpriv                  -0.012***         
                              (0.004)          
                                               
npub                         -0.011***         
                              (0.004)          
                                               
nelite                         0.006           
                              (0.005)          
                                               
hascollege                   -0.049***         
                              (0.006)          
                                               
popdensity2010              -0.0001***         
                             (0.00001)         
                                               
med_hhinc2016               0.00000***         
                             (0.00000)         
                                               
Constant                     0.251***          
                              (0.012)          
                                               
-----------------------------------------------
Observations                    741            
R2                             0.281           
Adjusted R2                    0.273           
Residual Std. Error      0.062 (df = 732)      
F Statistic           35.756*** (df = 8; 732)  
===============================================
Note:               *p<0.1; **p<0.05; ***p<0.01

nelite on kfr at different levels

htmtable.nelite <- stargazer(nelite.p1, nelite.p25, nelite.p50, nelite.p75, nelite.p100,
          type = "html",
          dep.var.labels = c("Bottom 1%", "25%", "50%", "75%", "Top 1%"),
          out = "nelitetable.html")

length of NULL cannot be changedlength of NULL cannot be changedlength of NULL cannot be changedlength of NULL cannot be changedlength of NULL cannot be changed

<table style="text-align:center"><tr><td colspan="6" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left"></td><td colspan="5"><em>Dependent variable:</em></td></tr>
<tr><td></td><td colspan="5" style="border-bottom: 1px solid black"></td></tr>
<tr><td style="text-align:left"></td><td>Bottom 1%</td><td>25%</td><td>50%</td><td>75%</td><td>Top 1%</td></tr>
<tr><td style="text-align:left"></td><td>(1)</td><td>(2)</td><td>(3)</td><td>(4)</td><td>(5)</td></tr>
<tr><td colspan="6" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left">nelite</td><td>-0.004</td><td>-0.002</td><td>-0.001</td><td>0.001</td><td>0.004</td></tr>
<tr><td style="text-align:left"></td><td>(0.003)</td><td>(0.003)</td><td>(0.002)</td><td>(0.002)</td><td>(0.003)</td></tr>
<tr><td style="text-align:left"></td><td></td><td></td><td></td><td></td><td></td></tr>
<tr><td style="text-align:left">popdensity2010</td><td>-0.0001<sup>***</sup></td><td>-0.0001<sup>***</sup></td><td>-0.00004<sup>***</sup></td><td>-0.00002<sup>***</sup></td><td>-0.00000</td></tr>
<tr><td style="text-align:left"></td><td>(0.00001)</td><td>(0.00001)</td><td>(0.00001)</td><td>(0.00001)</td><td>(0.00001)</td></tr>
<tr><td style="text-align:left"></td><td></td><td></td><td></td><td></td><td></td></tr>
<tr><td style="text-align:left">med_hhinc2016</td><td>0.00000<sup>***</sup></td><td>0.00000<sup>***</sup></td><td>0.00000<sup>***</sup></td><td>0.00000</td><td>-0.00000<sup>***</sup></td></tr>
<tr><td style="text-align:left"></td><td>(0.00000)</td><td>(0.00000)</td><td>(0.00000)</td><td>(0.00000)</td><td>(0.00000)</td></tr>
<tr><td style="text-align:left"></td><td></td><td></td><td></td><td></td><td></td></tr>
<tr><td style="text-align:left">Constant</td><td>0.221<sup>***</sup></td><td>0.350<sup>***</sup></td><td>0.471<sup>***</sup></td><td>0.585<sup>***</sup></td><td>0.770<sup>***</sup></td></tr>
<tr><td style="text-align:left"></td><td>(0.012)</td><td>(0.010)</td><td>(0.009)</td><td>(0.009)</td><td>(0.009)</td></tr>
<tr><td style="text-align:left"></td><td></td><td></td><td></td><td></td><td></td></tr>
<tr><td colspan="6" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left">Observations</td><td>741</td><td>741</td><td>741</td><td>741</td><td>741</td></tr>
<tr><td style="text-align:left">R<sup>2</sup></td><td>0.141</td><td>0.104</td><td>0.055</td><td>0.016</td><td>0.030</td></tr>
<tr><td style="text-align:left">Adjusted R<sup>2</sup></td><td>0.137</td><td>0.100</td><td>0.051</td><td>0.012</td><td>0.026</td></tr>
<tr><td style="text-align:left">Residual Std. Error (df = 737)</td><td>0.068</td><td>0.057</td><td>0.050</td><td>0.047</td><td>0.053</td></tr>
<tr><td style="text-align:left">F Statistic (df = 3; 737)</td><td>40.261<sup>***</sup></td><td>28.387<sup>***</sup></td><td>14.380<sup>***</sup></td><td>4.096<sup>***</sup></td><td>7.651<sup>***</sup></td></tr>
<tr><td colspan="6" style="border-bottom: 1px solid black"></td></tr><tr><td style="text-align:left"><em>Note:</em></td><td colspan="5" style="text-align:right"><sup>*</sup>p<0.1; <sup>**</sup>p<0.05; <sup>***</sup>p<0.01</td></tr>
</table>

Add a new chunk by clicking the Insert Chunk button on the toolbar or by pressing Ctrl+Alt+I.

When you save the notebook, an HTML file containing the code and output will be saved alongside it (click the Preview button or press Ctrl+Shift+K to preview the HTML file).

The preview shows you a rendered HTML copy of the contents of the editor. Consequently, unlike Knit, Preview does not run any R code chunks. Instead, the output of the chunk when it was last run in the editor is displayed.